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03 ; ABSTRACT 

Kepler observations of the star HD 187091 (KIC 8112039, hereafter KOI-54) revealed a remarkable light curve 

' exhibiting sharp periodic brightening events every 41.8 days with a superimposed set of oscillations forming a beating 

, pattern in phase with the brightenings. Spectroscopic observations revealed that this is a binary star with a highly 

i eccentric orbit, e=0.83. We are able to match the Kepler light curve and radial velocities with a nearly face-on (i=5.5 

' degree) binary star model in which the brightening events are caused by tidal distortion and irradiation of nearly 

t**"*- \ identical A stars during their close periastron passage. The two dominant oscillations in the light curve, responsible 

t-H . for the beating pattern, have frequencies that are the 91st and 90th harmonic of the orbital frequency. The power 

' spectrum of the light curve, after removing the binary star brightening component, reveals a large number of pulsations, 

, 30 of which have a signal-to-noise ratio > 7. Nearly all of these pulsations have frequencies that are either integer 

— < ■ multiples of the orbital frequency or are tidally-split multiples of the orbital frequency. This pattern of frequencies 

7-H ' unambiguously establishes the pulsations as resonances between the dynamic tides at periastron and the free oscillation 

r* \ modes of one or both of the stars. KOI-54 is only the 4th star to show such a phenomenon, and is by far the richest 

^ ■ in terms of excited modes. 

k> ■ Subject headings: binaries: close — binaries: spectroscopic — stars: individual (KIC 8112039, 

HD 187091, BD+43 3355, 2MASS J19461553+4356513) — stars: oscillations - 
03 ■ stars: variables: general 



1. INTRODUCTION 

In its search for transiting extrasolar planets, the Ke- 
pler Mission obtains high-precision time series photome- 
try of -156,000 stars (Koch 2010, Borucki 2010). This 
very large high-signal-to-noise sample, combined with 
the unprecedented near-continuous long-timescale cover- 
age, results in a high potential for discovery of astrophys- 
ically important stars. KOI-54 is one of those serendipi- 
tous discoveries. 

t Based in part on observations obtained at the W. M. Keck 
Observatory, which is operated by the University of California 
and the California Institute of Technology. 

To whom correspondence should be addressed: 
E-mail: wfw@sciences.sdsu.edu 



Cataloged as HD 187091 before being flagged as Ke- 
pler Object of Interest 54, this star was classified as a 
single bright (V=8.38) A star. The object did not seem 
unusual in any way prior to these Kepler observations 
and attracted little attention. This is in stark contrast 
to the astonishing Kepler light curve, which immediately 
offered a challenging puzzle begging for elucidation - 
see Fig. 1. 

In this paper, we present the Kepler observations and 
a large set of precise radial velocities that led to our bi- 
nary star model for KOI-54. We investigate the rich set 
of oscillations that are present, and conclude that these 
are tidally-driven pulsations. We also explore the evo- 
lutionary status of these stars and find a self-consistent 
scenario for the binary system. We finish with some re- 
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marks on the puzzles that still remain, and discuss in 
some detail the systematic issues that could potentially 
affect our interpretation, but find these issues to gener- 
ally have no significant effect on our conclusions. 

2. OBSERVATIONS 

2.1. Overview of the Light Curve 

KOI-54 has been observed by Kepler almost continu- 
ously from 2009 May 02 during Commissioning observa- 
tions (Q0) though 2010 Mar 22 "Quarter 4" (Q4). The 
light curve exhibits two remarkable features: a periodic 
brightening spike of ~0.7% occurring every 42 days and 
a ^0.1% "beat pattern" of pulsations in phase with the 
brightening events. The beat pattern arises from the in- 
terference of two pulsations with periods near 11 hours. 
The pulsations continue during the brightening and re- 
peat from event to event - they are strictly in phase with 
the brightening. 

An initial hypothesis that the repeating brightening 
events might follow from near-field microlensing (Sahu 
and Gilliland 2003) was pursued. Such a large event 
would result from a black hole of several solar masses 
orbiting and transiting an A star. Motivated by this 
hypothesis, a 1958 s duration SWIFT XRT observation 
was made on 2010 Apr 25, during one of the brightening 
events. The source was not detected in the 0.3-10 keV 
bandpass, but we estimate the 5-sigma upper limit on 
the X-ray flux to be F < 1.23 x 10~ 13 erg cm~ 2 s -1 . 
Using the HIPPARCOS parallax of 3.14 milliarcsec, the 
upper limit to the luminosity is L < 1.5 x 10 30 erg s -1 . 
Given the null detection, we make no further discussion 
of these observations. 

Radial velocity observations quickly ruled out the tran- 
siting black hole hypothesis because no large rapidly 
changing reflex motion of the A star was observed. Also, 
no ready explanation of the oscillations existed in this 
scenario, and when examined in detail the 42-day bright- 
enings deviate from the expected shape of microlensing 
events. Instead, as we show below, KOI-54 is a highly 
eccentric binary star system, and the brightening events 
are caused by tidal distortion and irradiation of the two 
stars during their close periastron passage. The two dom- 
inant pulsations producing the beating seen in the light 
curve have frequencies that are exactly 90 and 91 times 
the orbital frequency, hence result from tidally driven 
pulsation modes. 

2.2. The Kepler Photometry and Calibration 

The Kepler CCDs are read out every 6.54 s (6.02 s live- 
time) and co-added on-board. In Short Cadence (SC) 
mode the signal is combined to achieve approximately 
58.8 s sampling cadence, and in Long Cadence (LC) 
mode the data are binned to 29.424 min cadence. KOI-54 
is heavily saturated in these 6 s exposures, and blooming 
affects approximately 30 pixels of the 74-pixel aperture 
(approximately 7 pixels diameter plus two columns of 
25 pixels). However, the electrons are not lost, and be- 
cause of Kepler's stability, superb relati ve photometry is 
achievable — see IGilliland et all (poll for a discussion. 
For more details on th e design and perfor mance of the 
Kepler photometer see iKoch et al. I ( 20101) and Jenkins 
et al. (2010a,b). 

Each Kepler pixel spans ~4 arcseconds so the very 
large photometric aperture includes a few very faint 



(by comparison) background stars. The brightest star 
(KIC 8112007) has Kp=17.6 (4850 times fainter than 
KOI-54) , and ground-based images with 1 arcsecond see- 
ing show no star brighter than 18th magnitude within 5 
arcseconds. We conclude that background starlight con- 
taminates the light curve by less than 250 ppm (0.025%), 
and thus variations in the light curve are intrinsic to KOI- 
54. 

In addition to gaps due to the scheduled quarterly 
spacecraft rolls, safe-mode events and other spacecraft 
anomalies are present; cosmic rays and other noise 
sources also contaminate the signal. A simple automated 
sigma-threshold rejection method could not be used be- 
cause of the presence of the complex oscillation signal, 
so the light curve was carefully examined and obvious 
outliers were omitted by hand: 158 points were removed, 
leaving a total of 14,277 observations. This corresponds 
to a duty cycle of an astounding 92.5% over the span of 
10.6 months. 

The observations span 321.7 d, nearly eight complete 
41.8 d orbital cycles. For each spacecraft roll, the tar- 
get is on a different CCD with pixels of different sen- 
sitivity, and thus jumps in the light curve occur from 
Quarter to Quarter. Mean fluxes over the 4 Quarters 
span 1.01 - 1.07 xlO 10 e~/cadence. To correct for these 
changes, and the more troublesome safemode-related, 
pointing-adjustment, and other medium-timescale dis- 
continuities, the "RAW" pipeline processed light curved 
were detrended in the following way. First, the bright- 
ening events were masked in orbital phase from 0.9 to 
1.1 (i.e., 20% of the light curve). Then using gaps in the 
time series to define sections, the sections were fit with a 
low-order polynomial (typically cubic, though linear or a 
constant was used if appropriate). The polynomial was 
then subtracted from the time series, and the remain- 
der divided by the mean of the polynomial. A value of 
+1.00 was then added to this quotient giving a relative 
flux with a mean of unity. The pipeline uncertainty es- 
timates were boosted by a factor of 6, based on the rms 
scatter of the residuals of an early model fit. This scal- 
ing factor was chosen for simplicity, though each Quarter 
should have its own scaling (e.g. values of 5.0, 5.5 and 
7.3 for three different Quarters were found). The scat- 
ter in the residuals was significantly in excess of what 
was expected just from the error bars, but in hindsight 
this was due to using too few sinusoidal components in 
our initial modeling (see below). Slightly overestimating 
the uncertainties has no adverse affect on the estimation 
of the stellar and orbital parameters, and helps account 
for systematic noise (e.g. photometric trends after a safe 
mode or other anomaly) that is not included in the sta- 
tistical uncertainties. 

One month of Short Cadence data was obtained dur- 
ing Quarter 3: for t h e cha racteristics of Kepler SC data, 
see IGilliland et al. I (|2010D . These observations contain 
43,990 points and span 30.035 days. Examination of 
these SC data, and their power spectrum at higher fre- 
quencies, did not reveal any features not already well re- 
solved in the LC data, therefore only the Long Cadence 
data are used for our analyses. 

3 Data are available at the Multimission Archive 
(MAS T) at the Space Telescope Science Institute: 
http : //archive . stsci . edu/kepler/data_search/search.php 
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Finally, the relative timing precision is good to 0.5 s, 
but the absolute timing is uncertain by as much as 
6.5 s, and this systematic uncertainty should be added 
to the statistical uncertainty in the epoch of perias- 
tron T n in Table 2 — see Keple r Data Release 8 Notes 
(jMachalek fe Christiansenll2010h . 

2.3. Spectroscopy and Radial Velocities 

A variety of telescopes and instruments were used to 
provide moderate to high resolution spectroscopy for ra- 
dial velocity (RV) and spectral modeling. The spectra 
show double sets of absorption lines revealing the binary 
star nature of the system, and in Fig. 2 we show a section 
of a Keck HIRES spectrum where pairs of lines are clearly 
exhibited. The stars are of similar spectral type and lu- 
minosity. To determine the stellar parameters we carried 
out an analysis of 11 relatively clean Fe I and Fe II lines 
using the LTE spectral synthesis code MOOG (jSnedenl 
1973), modified for binary star analysis. Using an input 
linelist and two separate model atmospheres, the code 
computes a synthetic spectrum for each stellar compo- 
nent. Given a velocity separation and luminosity ratio, 
the code then overplots the resulting spectra onto the ob- 
served spectrum. We then varied the stellar parameters 
(temperature, gravity, metallicity and microturbulent ve- 
locity) until the best overall match to our selected iron 
lines was found. Using an HET HRS spectrum taken 
very near periastron (when the stars were cleanly sepa- 
rated in Doppler velocity by 25.5 km s _1 ), we obtained 
the following information about the stars. The stellar 
temperatures Ti and T 2 , are 8500 and 8800 K ±200 K, 
with log g of 3.8 and 4.1 ±0.2, and a luminosity ratio 
L2/L1 of 1.22 ±0.04 for the wavelength range covered 
in the spectral analysis, ^4500-6500 A. Both stars are 
~ 2-3 times more metal- rich than solar with [Fe/H] = 
0.4 ±0.2. Thus the stars are both A-type near-main se- 
quence star^aiidpotentially pulsating S Scuti variables 
— see lAerts et al. I ()2010[ ) for a thorough discussion of S 
Scuti and other pulsating stars. 

Estimates of V TO t sin i depended on the spectra and 
method used, but all indicate a low projected velocity 
ranging from < 5 to 10 km s _1 . Given the difficulty in 
measuring V ro t sini we adopt a value of 7.5 ±4.5 km s _1 . 
When corrected for the very low inclination of the sys- 
tem (discussed in §3) the true V rot can be < 50-100 
km s _1 , and this is still relatively low for an A-star. This 
low rotation rate, combined with the metal-rich abun- 
dance, suggests these are chemically peculiar Am or pos- 
sibly Ap stars. An inspection of the Nd III 6145A and 
Prill 6160A lines show that they are not strong, as they 
usually are in roAp stars. 

A total of 51 pairs of radial velocities were obtained 
in 2010 from six different telescopes±spectrographs, and 
these are listed in Table 1 . To minimize any potential sys- 
tematic offset between the velocities acquired with differ- 
ent instruments, all RV observations were calibrated us- 
ing the same procedure a nd all velocities measur ed using 
the TODCOR technique (jZucker fe Mazehlll994[ ). Every 
velocity measurement is referenced to a common radial 
velocity standard, HD 182488, that was observed every 
night, following standard practice for Kepler follow-up 
observations of targets of interest. Although not simul- 
taneous with the Kepler photometry presented in this 



paper, these observations provided the key to under- 
standing the nature of KOI-54: the stars are on a highly 
eccentric orbit, e=0.83, with periastron passage at the 
times of the brightenings. Thus the mutual interaction 
of the stars when closest together produces the bright- 
ening events seen in the Kepler photometry, and this is 
discussed in detail in §3 and §5.1. Fig. 3 shows the radial 
velocities and Keplerian fits, and Table 2 includes the 
orbital elements from the radial velocity-only fit. The 
ratio of K- velocities gives a mass ratio of 1.034, again 
confirming the similarity of the two stars. Note: The 
TODCOR methodology we employed for measuring the 
radial velocities does not provide uncertainty estimates 
on the radial velocities, so the mean of the residuals from 
the radial velocity fit to the observations was assigned as 
the uncertainty to all the velocities: 0.31 km s _1 . Given 
the relatively good fit to the RV data over most of the 
orbit, this approximation is justifiable. Concerns over 
systematic errors, as seen in the residuals of the fit, are 
discussed in §5.1.3. 

3. MODELING 

We employ a modified version of the E LC modeling 
code of Orosz (lOrosz fe Hauschi ldt 2000) to simultane- 
ously fit both the photometry and radial velocities. The 
non-spherical stars are covered with a fine grid of tiles, 
and for each time the intensity and velocity of the tiles 
are summed to produce the light curve and radial ve- 
locities. Gravitational distortions are modeled assum- 
ing a standard Roche potential, including the rotation of 
the stars themselves; the potentials are recomputed at 
each orbital ph ase because of the e l liptica l orbit — see 
the appendix of lOrosz fe Hauschi ldt ( 200 0]) for a descrip- 
tion o f the pote ntial that is ba s ed on lAvni fe Bahcalll 
(|1975h . and see iSepinksv et al. I ((20071 ) for a thorough 
discussion of equipotential surfaces in nonsynchronous 
eccentric binaries. Gravity darkening is included, using 
an ex ponent 0=0 . 25 appropriate for e arly- type radiative 
stars (jvon Zeipell 119241: iClaretl l2000h . We used tabu- 
lated spherical NextGen/PHOENIX model intensities 
(Hauschildt et al. 1997). The model flux is then inte- 
grated over the Kepler spectral response (approximately 
4250-8950 A, peaking at 5890 A with a mean wave- 
length of 6400 A — see lVan Cleve fc Caldwelll (|2009h and 
iKoch et al. I ((2010)). Irradiation of the stars is h andle d 
following the standard prescription of IWilsorj (1990): 
The light from each star is the sum of (1) the intrinsic 
intensities of each tile (modified for the local gravity, and 
if blackbodies are used, the limb darkening); and (2) an 
irradiation "reflection" component on the inward-facing 
hemispheres. The irradiation modifies the local temper- 
ature in the following way: T' 4 = T 4 x [l + Aboi^f^] 
where T and F are the temperature and bolometric flux 
of the star, and Fi rr is the incident bolometric flux from 
the companion star. Ab i is the bolometric albedo (not 
to be confused with the Bond albedo) and is the ratio of 
re-radiated to incident energy. A radiative atmosphere 
has Ab i = l (implying local en ergy conservation), a nd we 
hold A bo i=l in our models. iKallrath fe Milond ((19991 ) 
give an excellent description of Wilson's method and we 
refer the reader to that source for more details. At pe- 
riastron, the maximum change in temperature over the 
surface of the stars, including gravity darkening and ir- 
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radiation, is 88 K (=1%) for Star 1 and 61 K (0.7%) for 
Star 2. At apastron, the difference is only 0.5 and 0.4 K. 

The free parameters in the model are: the stellar 
masses Mi and M 2 , radii Ri and R 2 , temperatures Ti 
and T 2 , and 6 orbital parameters: inclination i, orbital 
period P, epoch of periastron T p , argument of perias- 
tron ui, eccentricity e, and systemic velocity 7 (held the 
same for both stars). The temperature is essentially un- 
constrained by the single-color broadband photometry, 
but the temperature ratio is weakly constrained by the 
ratios of radii and luminosities. While the models had 
freedom to vary the temperatures, the solution remained 
at the input spectroscopic temperatures. The individual 
radial velocity measurements are fit, not just the K\ and 
K 2 radial velocity amplitudes from the RV-only fit (i.e. 
we do not adopt the RV-only solution). The rotation of 
the stars, usually defined via the ratio fi = £l rot /£l or b, 
where Sl r ot is the stellar rotation angular frequency and 
£l or b is the orbital angular frequency, were also treated 
as free parameters. But for KOI-54 we defined these as 
£1 = Q rot /Q ps where f2 ps is the "pseudosynchronous" ro- 
tation frequency and depends only on the eccentricity 
and orbital period as given by eqn. (42) in Hut (1981): 

1 + 15 p 2 , 45 4 , _5_ e 6 

For an elliptical orbit true synchronous rotation is im- 
possible, but there is a pseudosynchronous spin such that 
over the course of an orbit there is no net torque on the 
star's rotation, and so the spin will not evolve: O = 
(Hut 1981). Note that Q ps is ~20% slower than the or- 
bital frequency at periastron, i.e., the spin is slower than 
what is necessary to keep the star tidally locked at peri- 
astron passage. 

In addition, other "observed parameters" are used to 
constrain the model: V rot smi, log g, and the ratio of 
luminosities L2/L1. These are not fixed in ELC, rather, 
the models arc steered toward them via a \ 2 penalty. 
A genetic algorithm and Markov chain Monte Carlo are 
used to find the best-fit models (in a \ 2 sense) and con- 
fidence intervals. The best-fit model has a reduced \ 2 
less than 1, indicating the uncertainties in the light curve 
were overestimated in the data calibration; but we deter- 
mined the l-<7 parameter intervals in the standard way by 
marginalizing over all other parameters and determining 
the interval bounded by Xmin + 1> we did not decrease 
this interval to account for the reduced \ 2 being less than 
1. 

3.1. ELC sinus 

The ELC code models the binary star light curve and 
radial velocity, but does not model the oscillations. At- 
tempts to isolate the binary light curve from the pul- 
sation light curve proved inadequate for the high preci- 
sion Kepler observations, so we added the following func- 
tionality into the ELC code, now dubbed "ELCsinus" . 
First, a trial binary star light curve is subtracted from 
the observed light curve. The mean is then subtracted 
and the Fourier transform is computed. The 15 largest 
peaks in the power spectrum are found (omitting side- 
lobes due to leakage of the 2 dominant peaks) , and a sum 
of sines and cosines is made using the 15 measured ampli- 
tudes and phases. This 15-sine pulsation model is then 



added to the binary star model to create the light curve 
model. In this ELCsinus model, the pulsations modulate 
the average flux and are not scaled by the binary star 
light curve; thus they are not boosted during the bright- 
ening events. The radial velocities are unchanged. The 
ELCsinus model thus consists of a light curve, a radial 
velocity curve for each star, and several derived param- 
eters: L2/L1, and the log g and V rot smi for each star. 
All of these are used to compute the goodness-of-fit \ 2 
statistic and parameter uncertainty estimates. 

As shown in Fig. 4, we obtain an excellent fit to the 
light curve, matching the amplitude, shape, and phase of 
the brightening events. The brightening is due to a com- 
bination of mutual irradiation and tidal distortion during 
periastron passage of the pair of A stars on their highly 
eccentric, nearly face-on orbit. The system parameters of 
KOI-54 are reported in Table 2. The uncertainties listed 
are the formal uncertainties and should be treated as 
lower limits, as there are many potential sources of sys- 
tematic errors as discussed in §5.5. As an example, the 
mass and radius of Star 1 from the ELCSINUS fit yields a 
log g that is 1.6 <r larger than the spectroscopically mea- 
sured log g (even though any deviation from the observed 
log g incurs a \ 2 penalty). Although it is very difficult 
to assess the size of systematic uncertainties, a doubling 
of the formal uncertainties seems reasonable, especially 
for the inclination, masses and radii of the stars. 

Note that the star designated as Star 1 is the less mas- 
sive star. Also notice that the pulsations continue to 
occur during the brightening events and exhibit a fixed 
pattern: there is no discernable phase drift. The pul- 
sations are matched well with our simple modification 
to ELC, indicating that the pulsations are not altered 
by the irradiation and gravitational distortions at pe- 
riastron. Indeed, as discussed in the next section, the 
pulsations are driven by events at periastron passage. 

4. PERIASTRON-PUMPED PULSATIONS 

To investigate the pulsations, the best-fit binary star 
model was subtracted from the light curve. This leaves 
the stellar pulsations, random noise, systematic calibra- 
tion noise, and any slight mis-matches between the model 
and actual binary star light curve. The 15-sines used in 
ELCsinus are not removed, only the binary star com- 
ponent of the model is removed. 

The beat pattern seen in the light curve indicates that 
two closely spaced pulsations are present in the light 
curve. This is confirmed by the power spectrum of the 
pulsation-only light curve - see Fig. 5. We computed 
the power spectrum using the Period04 software pack- 
age (Lenz & Breger 2005) which uses an iterative least- 
squares fit to all detected sinusoidal terms simultane- 
ously. Including the 2 dominant peaks in the power spec- 
trum, there are 13 clean peaks with amplitude greater 
than 12 ^mag, corresponding to a signal-to- noise ratio 
>20. We list the 30 strongest pulsations in Table 3. 
Frequencies less than 2 /xHz (longer than ^6 d) are not 
included in the Table as these long timescales suffer con- 
tamination from imperfect detrending, but we did in- 
clude 5 low-frequency terms when computing the power 
spectrum. The largest spike, Fl, is at a frequency of 
2.15286 d" 1 (24.9174 ^Hz) and the next largest, F2, at 
2.17680 d- 1 (25.1944 /iHz), or roughly 11.15 and 11.03 
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hours. These pulsations are perfectly sinusoidal; there 
are no harmonics and no evidence of any modulation in 
frequency (caused, for example, by Doppler shifts due to 
orbital motion). We cannot tell which star is pulsating, 
or if both are, and if so, which pulsations originate on 
which star. The pulsation phases are equal very near (but 
not exactly at) the time of periastron, thus the beating 
envelope amplitude is largest near periastron and lowest 
at orbital phase 0.5. These two pulsations are also the 
fastest in the light curve; there is no significant power at 
frequencies above F2. More precisely, between 3 d _1 (35 
^Hz) and the Nyquist frequency for the Long Cadence 
data at 24.5 d _1 (283 /iHz), there are no peaks above 3 
//mag. The one month of Short Cadence data was also 
examined, after simply omitting the single brightening 
event during this month, and its power spectrum also re- 
veals no signal greater than 3 /imag out to 720 d _1 (8.33 
mHz), showing a complete absence of any p-mode pul- 
sations. Thus neither star is a S Set star to exquisitely 
high precision, which given their early A spectral type, 
is somewhat surprising. 

The separation between the two largest peaks, in pe- 
riod, is 41.771 d, very close to the orbital period 41.805 d. 
Because the pulsations are not completely resolved in 
the power spectrum, their frequency difference is con- 
sistent with being identical to the orbital period. We 
assert that these are exactly the 90th and 91st multi- 
ples of the orbital frequency. This claim is confirmed by 
the other pulsations: 23 of the 30 largest pulsations are 
also very nearly exact integer multiples of the orbital fre- 
quency — see Table 3. Like the two dominant pulsations, 
these other pulsations are also very "pure" , with no mea- 
surable deviations from being perfectly sinusoidal. (The 
naming of the pulsations roughly corresponds to the rel- 
ative strengths of the pulsations, so F5 is the 5th largest 
pulsation. However, this is not an exact match because 
the amplitudes depend on the specific tapering method 
used, and as additional Quarters of data were added, 
some spikes swapped relative heights.) We searched for 
additional patterns in the power spectrum other than 
harmonics of the orbital frequency. This led to numer- 
ous detections of frequency spacings within multiplets at 
a value of Sf ~ 0.132 d -1 (^1.53 /iHz), e.g., between F6 
& F8, and F8 & F9 - see Table 3. Given this strong 
coupling of the orbital frequency with the pulsation fre- 
quencies, it is very likely that the pulsations are a result 
of a resonance between the dynamic tides and on e or 
more free low- frequency g-modes — see lAertd (|2007h and 
iWillems fc Aertsl ((2002) for a discussion of tidally-driven 
pulsations in binary stars. We will return to discuss this 
in §5.3. 

5. DISCUSSION, CONCERNS, AND SPECULATION 

5.1. Binary Star Phenomena 

KOI-54 is a remarkable binary star system. In the 
following section we further discuss various properties of 
the system, but first we mention that if the Hipparcos 
parallax is reliable for the binary, the distance is 318 ± 
71 pc. At apastron the stars are separated by 0.79 AU, 
which translates to 2.5 milliarcseconds - a potentially 
resolvable separation. 

5.1.1. Photometric determination of e audi 



The ELCsinus model produces an excellent fit to the 
photometry and tightly constrains the system parame- 
ters. To our knowledge, this is the first determination of 
stellar parameters for a non-eclipsing double-lined spec- 
troscopic binary based on the brightening during perias- 
tron passage. This is possible because of (i) the high pre- 
cision of the Kepler photometry and (ii) the strong sensi- 
tivity of the brightening to periastron passage, and hence 
to the orbital parameters. Remarkably, we have found 
that the photometry alone can determine the orbital ec- 
centricity as precisely as the radial velocities! This is a 
consequence of the dependence of the amplitude of the 
brightening to changes in eccentricity as illustrated in 
Fig. 6. In this set of simulations, closely matched to 
KOI-54, the only parameter that varies is the eccentric- 
ity. In particular, the lower panel shows how one can 
estimate the eccentricity to ~1% just from its ampli- 
tude, if all other parameters were known. The strong 
sensitivity to eccentricity is simply a consequence of the 
brightening being due to tidal forces and irradiation that 
scale as the separation of the stars cubed and squared, 
respectively, and the separation of the stars at perias- 
tron is linear in the eccentricity. (Note: we are assuming 
simple tidal distortion aligned along the line joining the 
center of masses; in the general asymmetric case the tidal 
force is much more sensitive to the separation, going as 
F oc r~ 7 — see Hut (1981) for a discussion.) 

In addition to the eccentricity, the Kepler photome- 
try is also able to constrain the orbital inclination, even 
though there are no eclipses and no double-humped ellip- 
soidal variations are present. But in fact the ellipsoidal 
variations are present - for such a highly eccentric or- 
bit the humps have shifted from the usual photometric 
quadrature phases, 0.25 and 0.75, to the phase of peri- 
astron (defined here as phase 1.0), and the humps have 
merged. The two humps are not equal in height because 
the orbital inclination is not exactly zero and the orbital 
ellipses are not aligned along our line of sight (the ar- 
gument of periastron is not ± 90 degrees); thus the two 
phases of maximum visible ellipsoidal distortion are not 
symmetric about periastron. In addition, the inclination 
and argument of periastron determine the orientation of 
the irradiated hemispheres of the stars to our line-of- 
sight, and this creates a small asymmetry and shift in the 
phase of the brightening. For the geometry of KOI-54, 
the larger the inclination, the more the peak would shift 
to earlier orbital phases, and the narrower the peak and 
the more asymmetric it would become (brighter post- 
peak than pre-peak). Such effects allow the inclination 
to be measured. 

5.1.2. Relative Contributions of Irradiation and Distortion 

Both tidal distortion and the irradiation/reflection 
contribute to the brightening, and Fig. 7 shows the sepa- 
rate contributions of each. For KOI-54 the non-sphericity 
is the larger contributor, with the reflection component 
roughly 75% as large. The dominance of the ellipsoidal 
effect is not a necessity, but rather a consequence of the 
specific binary system parameters, and could reverse in 
another system (i.e. reflection effect > ellipsoidal effect). 
As expected, the ellipsoidal component is more centrally 
peaked than the reflection component, as it is more sen- 
sitive to the separation of the stars. Note that the re- 
flection component is slightly asymmetric, being brighter 
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after periastron than before. This is a consequence of the 
inclination of the orbital plane and the orientation of the 
irradiated hemispheres with respect to the observer. 

5.1.3. Fit to the Radial Velocities 

While the ELCsiNUS model provides a good match 
to the Kepler photometry, the match to the radial ve- 
locities is not satisfactory at certain phases. The overall 
features are reproduced, and the rapid rise and fall of the 
velocities near periastron are well-fit, but the velocities 
at phases other than periastron are underestimated in a 
systematic way — see the residuals in Fig. 3. When we 
fit the radial velocities without the photometry, we ob- 
tain much better agreement: reduced x 2 =1.00 vs. 1.61 
for (102-8=94) degrees of freedom - the solid and dashed 
RV curves in Fig. 3 show the difference when using and 
not using the photometry. The maximum velocity differ- 
ence between the models is 0.3 km s _1 and occurs near 
orbital phase 0.1, agreeing with where the deviation be- 
tween the actual data and the ELCsiNUS fit are promi- 
nent. The fit is also poor near orbital phase 0.8, and 
it is unclear if between these two phases the fit would 
show the same systematic underestimate, since most of 
this phase range lacks observations. As a check, we fit 
the velocities using a well-tested RV-only code and ob- 
tained the same results. Also note that ELC computes 
an intensity-weighted net radial velocity, so that irradia- 
tion and tidal distortion are taken into account, e.g. cen- 
ter of light need not be the center of mass, especially at 
periastron. Thus the best-fit solution using the photom- 
etry is not the best-fit solution for the radial velocities 
alone. To confirm this, we used the best-fit radial veloc- 
ity parameters to compute a light curve, and the result 
was a very poor match to the photometry. The origin 
of the photometric versus radial velocity disagreement is 
not understood, and we can only speculate at this point. 
The poorest fit radial velocities are the ones with the 
lowest amplitudes, and these are the hardest to measure, 
as the lines from the two stars are most blended. This 
could account for most of the discrepancy. The stellar 
pulsations will affect the line profiles, thus adding jitter 
and broadening to the velocities which could induce a 
bias. The pulsations will also affect the measurement of 
V ro t sini, which could explain the wide range of V ro t smi 
we measure, from ^4 to 10 km s _1 . Finally, we have 
assumed that the stellar rotation axes are both aligned 
with the orbital axis. Although not expected to be mis- 
aligned (because the strong tidal forces will tend to align 
the spin axes, and the measured K-otsini are the same 
for both stars), a misalignment could induce significant 
jitter into the radial velocities and bias in the V rot sin i 
estimates. 

5.1.4. Tidal evolution 

With its eccentricity of 0.83, KOI-54 is very near the 
upper edge of the envelope of the distribution of eccen- 
tricities for its orbital period, e=0.89 following Mazch 
(2008). The "periastron period" , as introduced by Mazch 
(2008), is the period of a circular binary having the same 
semi-major axis as the periastron distance, and for KOI- 
54 Pperi-dist = 2.8 d. This short timescale suggests that 
the orbit should have experienced enough tidal evolution 
to synchronize over its main sequence lifetime and pos- 
sibly reduce its eccentricity from an even larger value. 



However, in his seminal work on the weak friction model 
for tidal interaction, Hut (1981, 1982) has shown that 
the evolution of eccentricity can be complex, and de- 
pends strongly on the equilibrium ratio of orbital and 
rotational angular momentum, a. In particular, for high 
values of a, the evolution of eccentricity need not be 
monotonically decreasing, as the spin and orbital angular 
momentum can be exchanged back and forth. Further- 
more, evolution of the stellar interior and its changing 
moment of inertia can alter simple orbital evolution (see 
Zahn 2008 for a discussion) . While there is over a factor 
of 10 uncertainty in a due to dependence on the un- 
known internal mass distribution (parameterized by the 
radius of gyration) and the poorly known stellar rotation 
period, the ratio of orbit-to-spin angular momentum is 
very large: ~ hundreds to thousands. This is well over 
a cr it=35.447, and tells us that the evolution of the orbit 
can be quite complex, with periods of increasing as well 
as decreasing eccentricity (Hut 1982). 

Tidal forces can synchronize the stellar spin on a 
timescale 2-3 orders of magnitude faster than circulariza- 
tion (e.g. Zahn 2008), so despite the large eccentricity, 
we can expect the stellar rotation to have experienced 
significant evolution towards synchronicity with the or- 
bit. The ratio of pseudosynchronous frequency to the 
orbital frequency is a function of the eccentricity only 
(Hut 1981) (see eq. 1), so we can compare this with the 
observations. The pseudosynchronous rotation frequency 
is fi ps =2.48 d _1 , translating to a spin period of 2.53 d, 
Vrot = 43-47 km s , and a projected V ro t sini=4.1-4.5 
km s . This value roughly agrees with the lower range 
of our observed V rot smi, and so a pseudosynchronous 
spin is roughly consistent with the observed e, P, and 
Vrot sinz. 

In our ELCsiNUS model we did let the stellar spins 
be free parameters; these are constrained mainly by the 
observed V ro t sin i but there is a weak dependence on the 
light curve because the oblateness of the stars depends 
on their rotation, and this affects the tidal distortion and 
reflection. In principle the stellar rotation can be refined 
through the light curve fit, and this could be used to 
check for pseudosynchronous rotation. But in practice 
this constraint is too weak to be useful. The models 
do have a marginal preference for Star 1 to have a ro- 
tation that is 3.5 times faster than pseudosynchronous, 
though the uncertainty admits the pseudosynchronous 
rate (there is a x 2 minima at 17 1 ~ 1), and such a rota- 
tion rate yields a V ro t sin i that is 2a high compared to 
the observation. If Star 1 is indeed spinning faster than 
pseudosynchronous, it is unclear if the star is still spin- 
ning down but has not yet reached pseudosynchronism, 
or if the star has exchanged angular momentum with the 
orbit and is being spun up from the pseudosynchronous 
rate. The observed slow rotation periods compared to 
normal A stars implies that both stars have experienced 
significant spin-down since birth; perhaps a spin-to-orbit 
exchange of angular momentum has prolonged or even 
increased the eccentricity of the orbit. However, Am and 
Ap stars may be born with a very slow spin, so this hypo- 
thetical exchange of spin and orbit angular momentum 
cannot be inferred. 

5.2. Exploratory Stellar Evolution Modeling 
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We have begun a search of parameter space for coeval 
stellar evolution models that fit the stellar constraints as 
given in Table 2. We are using an updated version of the 
Iben (f 963, 1965) stellar evolution code that includes the 
OPAL opacities (Iglesias fe R o gers 1 996) supplemented 
by the lAlexander fe Ferguson! ( fl994j ) and D. Alexander 
(1995, private communication) low-temperature opaci- 
ties, and the Grevesse & Noels (1993) solar element 
abundance mixture. As a starting point, we evolved a 
Mi=2.33 M© and a M 2 =2.40 M Q model with solar abun- 
dances: X=0.68 (H mass fraction), Y=0.28 (He mass 
fraction), and Z—O.OA (mass fraction of elements heavier 
than H and He). We find that the constraints in Ta- 
ble 2 are difficult to match with solar-metallicity Z=0.02 
models, especially the combination of temperature and 
radius: the radii are too small for the temperatures. A 
metal abundance higher than Z=0.04 would allow a wider 
range of parameters of the evolution models to match 
the (Mi,Ri,Ti)-(M2,R2,T 2 ) pair, and indeed the spec- 
troscopy does indicate enhanced metallicity. However, 
for this initial investigation we did not further pursue 
the metallicity as this would also require varying the Y 
abundance (beyond the scope of this investigation) and 
also that these stars may be metallic A stars (Am) so 
the surface abundances from spectroscopy may not be 
representative of the entire star. 

In selecting models, we emphasized the constraints 
on the ratios of temperature, luminosity, and radius 
over their actual values as the ratios are better deter- 
mined observationally. We discounted the spectroscopic 
log g measurements that suffer from large uncertainties. 
To be specific, we sought to match the ratios in tem- 
perature T 2 /Ti=1.035 ±0.034, luminosity L 2 /L x =1.22 
±0.04, and in radius R 2 /Ri=1.06 ±0.02. As expected, 
a family of solutions were found; for example a pair of 
stars with temperatures Ti, T 2 of 8498 K, 8667 K satis- 
fies our matching criteria, as does another pair at 8754 K, 
8999 K. Looking more closely, in the cooler pair, with an 
age of 0.354 Gyr, the temperature, luminosity, and radius 
ratios (1.020, 1.26, and 1.08) are close to the measured 
values, but the absolute radii (R 2 =2.52 and Ri=2.35 
Rq) are somewhat larger than our estimated values of 
2.33 and 2.20 R©. In the hotter pair, with an age of 
0.263 Gyr, the temperature, luminosity, and radius ratios 
(1.028, 1.25, and 1.06) are again close to the measured 
values, and the absolute radii (2.28 and 2.15 R©) are 
closer to the reported values, but the temperatures are 
at the high end of the uncertainty range. From this ex- 
ploratory investigation, the important conclusion we can 
draw is that the (Mi,Ri,Ti)-(M 2 ,R 2 ,T 2 ) set we measure 
is consistent with stellar evolution models of identical age 
and metallicity stars. 

5.3. Tidal Pulsation Modeling 

By way of brief introduction, linear perturbation the- 
ory can be used to determine the pulsation frequencies 
in stars, given their internal structure. Stellar pulsations 
are generally of two classes, the p-modes or g-modes de- 
pending upon whether pressure or gravity is the local 
restoring force. Being acoustic waves, p-mode character- 
istics are dependent on the local sound speed, and are 
important in the envelope and outer portion of the star. 
Their motion is primarily vertical and they have periods 
of order hours and less. The g-modes are set up by the 



competition between gravity and buoyancy deep within 
the star. Their wave motion is primarily horizontal and 
they have periods of order one day (for non-degenerate 
stars). The g-modes cannot extend through convective 
zones as the convective instability means there is no lo- 
cal restoring force. The modes are classified according 
to three spherical harmonic integers: the number of ra- 
dial nodes (n), the angular degree (£), and azimuthal or- 
der (to). ("Nodes" are the locations where the standing 
wave pattern, caused by the interference of the pulsa- 
tions, has zero amplitude.) The degree £ gives the num- 
ber of nodal circles on the sphere, and to is the number 
of nodal lines that cross the equator (or equivalently, go 
through the poles). Radial modes do not change the 
spherical symmetry of the star during the pulsation cy- 
cle and are characterized by £ = to = 0. The frequencies 
of non-radial modes of degree I are split by rotation or 
magnetic fields. Rotation lifts the degeneracy in m in 
such a way that each mode of degree £ results in 2£ + 1 
multiplet components in the frequency spectrum, with 
to = • • • , 0, • • • , +£. The classical k (opacity valv- 
ing) mechanism caused by He ionization in the stellar 
envelopes is predicted to drive p-mode pulsations in A 
to early F-type main sequence stars, while g-mode pulsa- 
tions are not expected in these ty pe stars. For a complete 
reference see lAerts et al. I ()2010l ). 

The periodic tidal forces experienced by the stars at 
periastron passage are expected to induce stellar oscilla- 
tions, including g-modes which otherwise would not be 
present. To explore such tidal excitation, we assume one 
of the stars (Mi) to rotate rigidly and the other star (M 2 ) 
to be well described as a point source. Whenever the ra- 
tio of the external tidal force to the self gravity at the 
star's equator is small, i.e. the dimensionless parameter 
£t = (i?i/a) 3 M 2 /Mi <C 1, the tide-generating poten- 
tial in the eccentric orbit can be expanded as a series 
in terms of spherical harmonics and of the companion's 
mean motion. We follow the method outlined in detail in 
iWillems fe Aertsl ((2002), who have shown that the spher- 
ical harmonic with degree £=2 is by far the dominant 
term in the expansion of the tide-generating potential, 
and only components with m = —2, 0, 2 occur. The val- 
ues of these components depend on the orbital inclination 
and eccentricity, as well as on the mean and true anomaly 
of the companion in its relative orbit. Whenever the forc- 
ing frequency, ctt = fc^orb + niQ ro t, matches the star's 
free gravity oscillation modes, the tidal action exerted by 
the companion gets enhanced and the particular oscilla- 
tion mode gets resonantly excited with the forcing fre- 
quency of the dynamic tide (Smeyers et al. 1998). Thus, 
we calculated the 1= 2, m=0 g-mode eigensolutions for 
each of these stars using the procedure and nonadiabatic 
non-radial pulsation code of Pesnell (1990) as reported in 
Guzik et al. (2000). We find that the g-mode spectrum 
for stars of this type is very dense, with period spac- 
ings of order 0.03 d. These are damped oscillations that 
would die out without the periodic tidal "pumping" , but 
the damping times are very much longer than the orbital 
period (on the order of 1000 years), so once a mode is 
excited, it is certainly not expected to change amplitude 
during an orbital period - or the lifetime of the Kepler 
Mission. 

While still adhering to the constraints on the stars dis- 
cussed above, we tweaked the tidal model parameters to 
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find a star whose I — 2, m = spectrum has one mode 
which coincides exactly with Fl (i.e., 90 times the or- 
bital frequency). This model has a mass of 2.33 M Q , 
T c ff = 8663 K and R — 2.15 Rq and most importantly, 
its <?i4 mode of t = 2, m = has frequency Fl, proving 
that standard stellar models exist that can generate a res- 
onant tidal excitation. (The gu mode has n=14 radial 
nodes in the interior of the star, and the £ = 2, m = is a 
quadrupole that distorts the star in the correct directions 
in response to tidal forcing.) Even though we did not 
tune the predictions of the tidal excitation theory to fit 
for KOI-54 in this discovery paper, we e xpect a frequency 
patter n as in the example treated in iWillems fc Aertsl 
(2002), i.e., a pattern of resonances at spaced values of 
the rotational and orbital frequency as fc£7 or b + mfl to t. 
We applied this to all the frequencies detected in KOI- 
54 and their splitting values, which allowed us to deduce 
that f2 ro t ~ 0.133 d _1 (1.54 //Hz) from the occurrence of 
m = 1 and m — 2 tidal splittings. This translates into a 
mean rotation period of 7.5 d and an equatorial rotation 
velocity of some 15 km s _1 . If the orbital and rotational 
axes are aligned, as assumed in the tidal theory, this 
leads to a prediction of V ro t sini of 1.5 km s _1 . This ro- 
tation velocity is certainly low compared to the observed 
T4. tsini, being ~ 1.7a from our adopted value, but we 
remark that the V ro t sin i measurements were made in the 
standard way and did not take into account pulsational 
broadening of the lines. We therefore tentatively accept 
the size of the tidal splitting of the lines as compatible 
with the measured V ro tsini. Furthermore, it is unlikely 
that the rotation inside the star will be rigid, and so this 
value of ftrot should be interpreted as an average value 
over the depth of the star, effectively negating concerns 
that this average rotation velocity is too low compared 
to the observed surface V ro tsm.i. In summary, we have 
found a model consistent with the observations that can 
explain all but one of the 30 largest observed pulsations 
as either harmonics of the orbital frequency or m = 1 and 
rn = 2 tidal splittings of those harmonics. (F25 remains 
unaccounted for in this model.) Despite the very attrac- 
tive features of this model, it is still very exploratory in 
nature and we do not claim uniqueness of the solution. 
Future observations will lead to much better frequency 
resolution, allowing a much more refined stellar model. 
We point out that slight deviations of the observed fre- 
quencies of the splittings from their exact predicted fre- 
quencies is quite interesting, and probably related to the 
nonrigidity of the stellar rotation. When several years of 
Kepler observations are available, we have the exciting 
prospect of using the observed tidal splittings to allow 
for an asteroseismic determination of the internal rota- 
tion law of the star. 

5.4. Pulsations Puzzles 

From their spectroscopic temperatures (T e g ^8500, 
8800 K) both stars are too hot to be 7 Dor pulsators, 
and the absence of any periodicities shorter than 11 h in 
the Short Cadence light curve power spectrum confirms 
this. Why the stars do not exhibit any 8 Scuti-like pul- 
sations is not known, though it might be that the stars 
are just beyond the blue edge of the classical instability 
strip. Presumably no unstable modes (p or g) exist for 
the stars, and the stars are intrinsically non-pulsators. 
This supposition is supported by the exploratory stellar 



models discussed in the previous section. 

The pulsations that are observed are those modes for 
which their free g-mode eigenfrequency is close enough to 
a harmonic of the orbital period that they are resonantly 
driven to measurable amplitudes by the periodic tidal 
(and possibly thermal) forces felt at periastron. Thus 
while very many pulsation modes are possible, only those 
that are integer multiples of the orbital frequency and 
phase-locked with the binary orbit are present. 

Some challenges posed by the pulsations in KOI-54 in- 
clude: (1) Why are the 90th and 91st harmonics of the or- 
bital frequency the strongest pulsations? (2) Why should 
one particular mode be strongly excited while many ad- 
jacent ones are not, i.e., why is orbital harmonic 91 very 
strong but 92 completely absent? (3) What limits the 
highest frequency excited mode and why is there no sig- 
nificant power present at higher harmonics? (4) Why are 
some modes tidally split but not others? (5) If only one 
of the stars is pulsating, which one is it, and why isn't 
the other pulsating? 

Interestingly, if we accept the slow stellar rotation in- 
ferred from the tidal splitting of the pulsations, that 
would mean the star is rotating at only a third of the 
pseudosynchronous rate. This would seem unlikely in a 
simple tidal evolution scenario, but the very large value 
of orbit-to-spin angular momentum ratio means that ex- 
changes of angular momentum are possible and episodes 
of spin-up and spin down can take place. Or, as men- 
tioned earlier, the stars are possibly Ap stars that may 
have been born rotating very slowly. A factor of two 
improvement in the measured V ro t sin i would be very 
valuable; the precision is currently available, but the ac- 
curacy is not. 

Given that the stars are similar in mass, radius, and 
temperature, have identical ages and metallicities, and 
experience identical driving frequencies, a total lack of 
pulsations from one of the stars would be puzzling. And 
yet the matching of the pulsation frequencies with a sin- 
gle pulsating star, and the lack of any apparent double- 
set of peaks in the power spectrum, does suggest that one 
star could be responsible for all the pulsations. In princi- 
ple, high-spectral resolution, time-resolved spectroscopy 
can determine which of the stars is pulsating (or if both 
are). If the spectra are obtained near periastron, the 
component lines will be well separated and pulsations in 
one of the sets of lines immediately tells us which star is 
pulsating. However, the intensity amplitude of the pul- 
sations is very small, ~ 0.1%, making this rather chal- 
lenging. Velocity modulations causing line profile shape 
changes are much more likely to be observable, though 
again the expected variations will be small. Further com- 
plicating the matter is the need to cover at least one full 
oscillation cycle, thus requiring >12 hours, and perhaps a 
multi-longitude campaign. Fortunately, the star is bright 
(Kp=8.380) and the integrations can be long and still 
temporally resolve the pulsations. The reward for such 
efforts is high: given that the two stars are very similar, 
if only one star is pulsating, it means that the excitation 
mechanism must be very fine-tuned to drive the oscilla- 
tions in one star but not the other. Thus there is the 
potential for a precision investigation of internal stellar 
structure, given stars of identical age and metallicity, and 
nearly identical masses. 
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5.5. Additional Systematic Concerns 

While we are confident the periastron-pumped pulsat- 
ing binary star model is correct, given the complexity of 
both the physics and the data calibration, several sources 
of systematic error are potentially present. In this sec- 
tion, we describe several of these, with an assessment 
of their effects when possible. These are checks of the 
robustness of the overall solution, not detailed investiga- 
tions that could reveal small changes at the few a level. 
Towards that goal, we used a linear limb-darkened black- 
body approximation instead of stellar atmospheres, and 
locked the stellar rotation to the pseudosynchronous spin 
frequency. We begin with a look at the data calibration. 

The most significant of the systematic errors in the 
data arise from the difficulty in removing small and 
abrupt changes in the flux level due to changes in point- 
ing, cosmic ray events, safe modes, etc. While such 
breaks are usually visible in Kepler light curves, the pul- 
sations in KOI-54 hide these. The residuals of the ELC- 
SINUS fit exhibit small but significant meanderings and 
tilts on longish timescales (days to weeks) which are al- 
most certainly due to detrending problems. Fig. 4 shows 
the residuals of the data minus the ELCsinus fit scaled 
by a factor of 5 to be more visible. The residuals are not 
white noise, indicating correlations due to the detrend- 
ing, but also due to the fact that the model includes only 
15 sinusoids while at least 30 are present. The largest set 
of correlated residuals occurs at the brightening event 
near day BJD 2,455,105 (see middle left panel in Fig. 4). 
Because the determination of the eccentricity and incli- 
nation are so sensitive to the amplitude and shape of the 
brightenings, even small tilts or jumps could potentially 
affect these estimates. We were concerned about the in- 
clination in particular: Normally (e.g. for eclipsing bina- 
ries) a small change in inclination does not significantly 
affect the mass estimates for the stars, as the slope of the 
sine function near 90 degrees is flat. But for KOI-54 the 
inclination is small so changes in sini go linearly with i. 
Because the masses depend on sin -3 i a change from 5 de- 
grees to 6 degrees results in a decrease of 58% in mass! To 
see if detrending imperfections could affect the accuracy 
of our model parameters (as opposed to precision), we 
carried out two investigations. First, we masked out the 
region with the largest residuals and re-fit. Second, we 
used the power spectrum feature in ELCsinus to mea- 
sure, then remove, low-frequency power in the residuals, 
on timescales of 10 d and longer. While the % 2 of the 
models was significantly better (by construction), there 
was no significant change in the system parameters. 

The determination of the eccentricity depends on the 
relative flux increase of the brightening event. Thus any 
systematic error in normalization will affect the eccen- 
tricity. The light curve is in fact contaminated by back- 
ground starlight contained in the same aperture used for 
KOI-54, and the relative amplitude of the brightening is 
thus biased low, and hence the derived eccentricity is po- 
tentially biased low. To check the effect this may have, 
we subtracted 250 ppm from the light curves and re-fit. 
In addition, we included a small "third light" dilution 
factor into the model and fit for this parameter. In both 
cases, there was no significant change in the model pa- 
rameters. Furthermore, such a bias would not effect the 
radial-velocity only eccentricity estimate, which agrees 



remarkably well with the photometric estimate. 

The "wings" of the brightening event (i.e. the upward 
curvature of the light curve) extend all the way to apas- 
tron at phase 0.5; there is no phase where the light curve 
is truly flat. This is a potential problem in the cali- 
bration: we detrended the light curve assuming it was 
flat between phases 0.1 and 0.9. From our binary-star 
only light curve, the maximum difference in relative flux 
between phases 0.1 through 0.9 is 99 ppm, roughly 1.7 
times the adopted uncertainty per point. However the 
light curve is flat to within 25 ppm from phase 0.2 to 
0.8, so the neglect of curvature in the calibration is be- 
nign. 

In our model, we have neglected several physical ef- 
fects which we address below. Systematic errors in the 
model are potentially more serious than those of data cal- 
ibration, as we have 8 brightening events that can help 
average over any calibration errors. For example, as is 
standard practice, we have ignored any consideration of 
magnetic fields. However, for these Am or Ap stars, such 
neglect at periastron passage may not be entirely justi- 
fied and warrants further investigation. 

We have not included relativistic Dopper boosting of 
the light curves. For this nearly face-on binary the 
boostin effect will be very small, and a quick test showed 
that including the boosting produced a change in total 
(not reduced) % 2 of ~0.1, a completely negligible amount. 
Nevertheless, in general, future models should include 
this effect as a modestly higher orbital inclination would 
result in a measurable difference in the light curve, given 
Kepler's remarkable precision. 

In general, an irradiated star is hotter than an un- 
irradiated star, and so the light that impinges back onto 
the irradiating star is greater than if the irradiated star 
were isolated. Thus the reflection effect should be itera- 
tively com puted, and thi s is the standard method cham- 
pioned by IWilsonl ()1990[ ). However, for KOI-54 this is 
a very small effect as the total brightening is less than 
1%, and tests with iteration showed almost no differ- 
ence compared to a single calculation assuming a point 
source irradiator (which is exact for perfectly spherical 
stars). Given the large computational cost of full geom- 
etry tile- by-tile iteration, and the negligible benefit, we 
used the point-source irradiation approximation. Our 
models have a maximum non-sphericity of 0.7% at peri- 
astron (ratio of minimum to maximum radii: pole radius 
to Ll-direction radius), justifying the spherical approxi- 
mation for irradiation. And despite the relatively close 
periastron passage of 0.065 AU (~6.4 stellar radii), the 
stars only fill < 30% of their Roche lobe radii. 

The most significant of the limitations of our model is 
the use of the Roche approximation. While quite success- 
ful at matching the observations, there is a systematic 
bias inherent in the use of Roche potentials in that the 
tidal elongations are always pointing along the axis con- 
necting the centers of mass. In other words, the stars are 
treated as perfect, viscousless fluids that can instanta- 
neously adjust their shapes to the gravitational and spin 
potentialsjj But real stars will have tidal bulges that do 
not instantly re-adjust to the changing external poten- 
tial, and are also "dragged" by the rotation of the star. 

4 The Roche potential approximation is similar to the simplest 
equilibrium tide models in this regard. 
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Thus the elongation will not be aligned with the cen- 
ter of masses except when the stars are in synchronous 
rotation in a circular orbit. The irradiation/reflection 
would also be affected, simply because of the different 
geometry. Furthermore, the Roche potential treats the 
mass distribution as a point-mass at the center of the 
star; this is equivalent to an n=5 polytrope or a tidal 
Love number k2=0. Thus for a given external force, the 
Roche potential yields the least possible tidal deforma- 
tion. It then follows that the amplitude of the ellipsoidal 
variation is minimal in the Roche approximation. While 
it is true that the mass inside stars is highly centrally 
concentrated, a more realistic mass distribution would 
produce a larger tide for a given external force. So to 
match a given observed brightening amplitude, a weaker 
tidal force would be required. This could be achieved by 
either lower masses, or much more likely, a very small 
increase in separation at periastron due to a very small 
decrease in eccentricity. Hence the brightening in the 
Roche potential model is expected to be very slightly 
different than a case where a more realistic internal mass 
distribution is used. 

Although we have measured very precise system pa- 
rameters with our ELCsinus model, for the reasons dis- 
cussed above the accuracy is worse than the precision. 
This is especially true for the eccentricity and inclination, 
and therefore the stellar masses which are proportional 
to the cube of the inclination. On the other hand, the 
photometric-only eccentricity agrees well with the radial 
velocity-only eccentricity, indicating that the Roche ap- 
proximation is valid at the level of analysis presented in 
this discovery investigation of KOI-54. A better treat- 
ment of the system, well beyond the scope of this paper, 
would involve dynamical tide theory — particularly rel- 
evant given the observed tidally-driven pulsations. 

This deficiency of the Roche model leads to an inter- 
esting supposition. The angle between the instantaneous 
centers of the stars and the axis of the tidal bulge de- 
pends on the internal structure of the star and is related 
to the efficiency of tidal dissipation. Measurement of this 
lag angle can then be used to constrain the tidal quality 
factor Q. If one could determine the instantaneous ge- 
ometry of the system without the photometry, then the 
lag between the brightening and the true time of peri- 
astron passage could be measured. The geometry could 
in principle be determined though a great host of radial 
velocities, corrected for the distortions caused by the pul- 
sations. Another possible way is related to the pulsations 
themselves; they provide a clock that could be used to 
set a fiducial time against which to measure the time of 
brightening. Interestingly, the time when the two largest 
modes (Fl and F2) are in phase is not the time of pe- 
riastron as defined by the photometry, but 0.39±0.07 d 
(=9.4 h) earlier. This can actually be seen in the light 
curves themselves, as the asymmetric shoulders of the 
brightening and dip at the very pinnacle - see Figs. 1 
and 4. 

5.6. KOI-54 in Context 

While tidally-excited oscillations have been studied in 
quit e some detail from a theoretical point of view (e.g., 
see (jAerts et al. II2010T ) for a discussion and references), 
observational evidence of their existence is scarce. Prior 
to KOI-54, only 3 such systems were known to exist. The 



rarity of these stars is understandable if one realizes that 
the forcing frequencies from the dynamical tides must 
come very close to the free eigenmode frequencies of one 
of th e stars in order to resonantly excite the oscillation, 
e.g., IWillemsl (|2003f) . This naturally leads to the exci- 
tation of particular gravity modes whose frequencies are 
exact multiples of the orbital frequency. 

Two of the previously known cases, HD 177863 (De Cat 
et al. 2000) and HD 209295 (Handler et al. 2002) were dis- 
covered from ground-based photometry and spectroscopy 
and are single-line spectroscopic binaries. HD 177863 is a 
slowly pulsating B star with e = 0.60, P or b = 11.9 d with 
one detected gravity mode whose frequency is 10 times 
the orbital frequency. HD 209295 is a hybrid 7 Dor/ 8 Set 
star in a binary with e = 0.35, P or b = 3.1 d with five 
g-modes having frequencies which are exact multiples of 
the orbital frequency, besides a free 6 Set mode. A much 
more interesting case is the eclipsing double-lined spec- 
troscopic B-type binary HD 174884 discovered with the 
C0R0 T mission (Maceroni et al. 2009). This system has 
e=0.29, P or b = 3.66 d and has pulsations with exact mul- 
tiples of 2, 3, 4, 8, and 13 times the orbital frequency and 
very tightly constrained physical parameters. While it is 
unclear if the lowest-order multiples are due to imperfect 
prewhitening of the orbital curve (i.e., removal of the bi- 
nary star contribution to the light curve) one does not 
expect this for harmonics 8 and 13 as all lower-order har- 
monics should have been found as well. Moreover, these 
two higher frequencies correspond exactly to those of free 
gravity modes of radial order ~ 10, which points to tidal 
excitation. 

With its periodic tidal brightening and rich set of pul- 
sations (over 30 pulsations at either integer multiples of 
the orbital frequency or at tidally split multiples of the 
orbital frequency are present), KOI-54 now joins this elite 
set of eccentric binaries that exhibit tidally-driven pul- 
sations. 

6. SUMMARY 

Far from being an ordinary A star, the exquisite Ke- 
pler observations of KOI-54 have revealed the object to 
be a fascinating binary star system exhibiting a host of 
interesting phenomena. We have successfully matched 
the light curve and radial velocities with a model con- 
sisting of a pair of A stars on a highly eccentric (e=0.83) 
orbit, seen nearly face-on (i=5.5 degrees). As the stars 
closely pass each other at periastron, coming within ^6 
stellar radii, the stars tidally distort each other's shape 
and mutually irradiate and heat each other. The com- 
bination of the tidal ellipsoidal variation and the irradi- 
ation/reflection effect creates the periodic 0.7% bright- 
ening seen every 41.8 days. In addition, the close pe- 
riastron passage is responsible for exciting a rich set of 
stellar pulsations with at least 30 modes. The two largest 
pulsations, at the 90th and 91st harmonic of the orbital 
frequency, beat against each other producing the mod- 
ulation envelope seen in the light curve. The remaining 
pulsations are explained as additional harmonics of the 
orbital frequency or tidally-split harmonics of those fre- 
quencies, clearly establishing the tidally-driven origin of 
the pulsations. Using these frequencies we are able to 
deduce the average rotation period of the pulsating star, 
7.5 d. KOI-54 is the only case where the tidally-excited 
oscillations have allowed such a measurement. Future 
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Kepler observations will allow factors of several times 
higher precision frequency measurements, offering the 
potential to map the internal rotation profile of one of 
the KOI-54 stellar components. 

Despite our success in modeling the system, many puz- 
zles remain. For example, it is not known which of the 
stars is pulsating, or if both are. Given that the stars 
are similar in mass and radius, and identical in age and 
metallicity, it is not known why one star would pulsate 
and the other not. Finally, although the pulsations and 
brightening events appear quite dramatic in the light 
curve, the amplitude of the pulsations is in fact consid- 
erably less than 0.1%; the discovery of the nature of this 
remarkable star system was made possible by the extraor- 
dinary precision and duration of the Kepler photometry. 

Kepler was selected as the 10th mission of the Dis- 
covery Program. Funding for this mission is provided 
by NASA, Science Mission Directorate. The authors ac- 
knowledge support from the Kepler Participating Scien- 
tists Program via NASA grant NNX08AR14G. CA and 
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the California Institute of Technology, the University of 
California, and the National Aeronautics and Space Ad- 
ministration. The Observatory was made possible by the 
generous financial support of the W.M. Keck Foundation. 
We thank Mr. Cur Windmiller for general assistance and 
for a careful reading of this manuscript. We especially 
thank the many members of the Kepler Team whose hard 
work made these observation possible. Finally, we thank 
the anonymous referee for a thorough review of this pa- 
per. 
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Fig. 1.— 



upper: The detrended and normalized Kepler light curve of KOI-54. lower: A detailed view of a brightening event. 
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Fig. 2.— 



Keck HIRES spectrum showing the double-line absorption and revealing the binary nature of KOI-54. 
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Fig. 3. — upper: Observed radial velocities (RV) and fits: the dashed line is the fit to the RV data only, the solid line is the best fit 
model that simultaneously fits the RV and the light curve. The filled circles denote the RV curve of Star 1. middle: The residuals of the 
data minus the best fit. lower: The ELCsiNUS model light curve showing the predicted fluxes at the time of the RV observations. 
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Fig. 4. — Light curves showing individual brightening events with the observations plotted as red points, the ELCSINUS model fit in 
black, and the residuals (offset by +0.998 and scaled by a factor of 5) in green. The lower right-hand panel shows just the model. 
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Fig. 5. — Power (amplitude) spectrum of the pulsation-only light curve, with the largest pulsations labeled. The bottom panel shows 
the power spectrum after prewhitcning by removing the two dominant pulsation Fl and F2. 
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KOI-54 Radial Velocities 
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Note. — Velocities reported in 
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Uncertainties arc estimated to be 0.31 km s _1 . 





a Facility Code: 

TRES: Tillinghast Reflection Echcllc Spectrograph on the F.L. Whipple Observatory 1.5-m telescope. 
FIES: FIbrc-fcd Echcllc Spectrograph on the Nordic Optical Telescope. 

MCD: Tull spectrograph on the McDonald Observatory 2.7-m Harlan J. Smith Telescope. 
HIRES: HIRES spectrograph on the W.M. Keck Observatory Keck I telescope. 
LICK: Hamilton Echcllc Spectrograph on the Lick Observatory Shane 3-m Telescope. 
HET: HRS spectrograph on the McDonald Observatory Hobby-Ebcrly Telescope. 
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TABLE 2 
KOI-54 System Parameters 
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Orbital eccentricity, e 


0.8335 


0.0005 




Arg. periastron, lu 


36.70 


0.90 


degrees 


Orbital inclination, i 


5.50 


0.10 


degrees 


Semimajor axis, a 


0.3956 


0.008 


AU 


Star 1: Qi 


3.5 


2.3 




Star 2: f_ 2 


1.0 


0.9 




Star 1 mass: Mi 


2.33 


0.10 


M 


Star 2 mass: M 2 


2.39 


0.12 


M 


Star 1 radius: Ri 


2.20 


0.03 


R 


Star 2 radius: R 2 


2.33 


0.03 


Rq 


Note. — O is defined as 


the ratio of rotation frequency to the pscudosynchronous 


rotation frequency: O = Q rot /£l P s 
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TABLE 3 



Thirty Largest KOI-54 Pulsations 


ID 


frequency 


frequency 


amplitude 


f 1 f orbit 


nearest 




(d- 1 ) 


(/iHz) 


(jumag) 




harmonic 


Fl 


2.1529 


24.917 


297.7 


90.00 


90 


F2 


2.1768 


25.195 


229.4 


91.00 


91 


F3 


1.0525 


12.182 


97.2 


44.00 


44 


F4 


0.9568 


11.074 


82.9 


40.00 


40 


F5 


0.5363 


6.207 


82.9 


22.42 


— 


F6 


1.6405 


18.988 


49.3 


68.58 


— 


F7 


1.7222 


19.933 


30.2 


72.00 


72 


F8 


1.5087 


17.462 


17.3 


63.07 


63 


F9 


1.3773 


15.941 


15.9 


57.58 


— 


F10 


0.6697 


7.751 


14.6 


28.00 


28 


Fll 


1.2678 


14.673 


13.6 


53.00 


53 


F12 


1.1241 


13.011 


13.4 


46.99 


47 


F13 


0.9329 


10.798 


12.5 


39.00 


39 


F14 


1.4349 


16.608 


11.6 


59.99 


60 


F15 


0.8851 


10.244 


11.5 


37.00 


37 


F16 


1.6983 


19.656 


11.4 


71.00 


71 


F17 


0.6183 


7.156 


11.1 


25.85 


— 


F18 


1.8178 


21.039 


9.8 


75.99 


76 


F19 


0.8574 


9.924 


9.3 


35.84 




F20 


0.6458 


7.475 


9.1 


27.00 


27 


F21 


1.0284 


11.903 


8.4 


42.99 


43 


F22 


1.0765 


12.460 


8.3 


45.01 


45 


F23 


1.5092 


17.467 


8.1 


63.09 


63 


F24 


0.8610 


9.965 


6.9 


35.99 


36 


F25 


1.4452 


16.726 


6.8 


60.42 




F26 


1.2439 


14.397 


6.4 


52.00 


52 


F27 


1.0078 


11.664 


6.3 


42.13 




F28 


0.7894 


9.137 


5.9 


33.00 


33 


F29 


0.6937 


8.028 


5.8 


29.00 


29 


F30 


1.1483 


13.290 


5.7 


48.00 


48 



Note. — Uncertainty in frequencies is ~0.0001 d 1 or 0.001 /iHz. Formal uncertainty in amplitudes is 0.3 fim&g. Orbital frequency f orbit was 
found via least-squares fit to best match the harmonics: f or bit — 0.0239205 d _1 — 0.276857 fiHz. 



